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Verfahren zur Erzeugung einer Einhiillenden 



Die Erfindung betrifft ein Verfahren zur automat ischen Erzeu- 
gung einer Einhiillenden fur ein rechnerverf iigbares Konstruk- 
tionsmodell, das durch eine vorgegebene Vernetzung in Finite 
Elemente zerlegt ist. 

Unter einer Einhiillenden wird eine rechnerverf iigbare und na- 
herungsweise Beschreibung der Oberflache des vorgegebenen 
Konstruktionsmodells durch eine Menge von Flachenelementen 
verstanden. Die Einhiillende approximiert die Oberflache mit 
einer vorgegebenen Maximal abweichung und vernachlassigt Ld- 
cher, Aussparungen, Bohrungen, Taschen, Rillen und andere Un- 
regelmafcigkeiten, wenn deren Abmessungen kleiner als die vor- 
gegebene Maximal abweichung sind. 

Eine solche Einhiillende wird benotigt, urn physikalische Vor- 
gange an der Oberflache eines technischen Systems durch eine 
Finite-Elemente-Simulation nachzubilden. Hierfiir wird das 
technische System durch ein rechnerverfiigbares Konstruktions- 
modell nachgebildet. Beispielsweise ist das modellierte Sys- 
tem ein Motor oder eine Karosserie eines Kraf tf ahrzeugs , und 
die Gerauschemission von der Oberflache des Motors bzw. der 
Karosserie wird durch eine Akustik-Simulation vorhergesagt . 
in der Akustik-Simulation werden nur Vorgange auf der Ober- 
flache und der angrenzenden Umgebung beriicksichtigt , nicht a- 
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ber Vorgange im Inneren des Motors bzw. der Karosserie. Damxt 
die Rechenzeit der Akustik- Simulation moglichst gering ist, 
approximiert die Einhullende die Oberflache nur so genau wie 
fur die Simulation notig. 

Ein Software-Werkzeug fur Akustik- und andere Simulationen 
namens „Virtual .Lab" wird unter http://www.lmsintl.com, abge- 
fragt am 24. 10. 2003, beschrieben. Erwahnt wird die Erzeu- 
gung einer Einhullenden, dort „ acoustic surface meshing pro- 
cess" genannt, aus einer vorgegebenen Vernetzung. Aus einer 
Vernetzung eines beispielhaf ten Antriebstrangs mit 150.000 
Finiten Elementen wird eine Einhullende mit 5.000 Flachenele- 
menten erzeugt . Nicht offenbart wird, wie die Einhullende er- 
zeugt wird. 

In DE 10023377 C2 wird ein Verf ahren zur Finite-Elemente- 
Simulation beschrieben. Ein vorgegebenes zwei- oder dreidi- 
mensionales Simulationsgebiet wird mit quadratischen Flachen- 
elementen bzw. wiirf elf ormigen Volumenelementen uberdeckt . Er- 
mittelt wird, welche der Gitterzellen ganz, teilweise oder 
gar nicht im Simulationsgebiet liegen. Wie dies geschieht, 
wird nicht offengelegt. Im Simulationsgebiet werden anschlie- 
fcend approximierende Funktionen mit Hilfe von B-Splines ilber 
den Gitterpunkten definiert, und die Simulation wird mit Hil- 
fe dieser Funktionen durchgef uhrt . 

Aus DE 10056107 Al 1st ein Verfahren zur Ermittlung von Klap- 
pergerauschen an Kraft fahrzeugen bekannt . Ein Konstruktions- 
modell eines Fahrzeugs wird in Elemente und Gitter unter- 
teilt. Weiterhin wird der Def initionsbereich des Modells in 
Unter-Definitionsbereiche unterteJ.lt, und die Elemente und 
Gitter eines Unter-Def initionsbereichs werden jeweils separat 
untersucht . 

in DE 19919891 Al wird die Erzeugung eines virtuellen Windka- 
nals z. B. fur ein Kraf tf ahrzeug beschrieben. Ein vorgegebe- 
nes Konstruktionsmodell wird vernetzt (eine „Masche" wird er- 
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zeugt), beispielsweise ein Oberf lachenmodell mit 60.000 drei- 
eckigen Flachenelementen. Fur diese Vernetzung werden vor- 
zugsweise kommerzielle Werkzeuge verwendet . Bei Bedarf wird 
das Konstruktionsmodell verfeinert, z . B . in Bereichen mit 
Verwirbelungen . 

Der Erfindung liegt die Aufgabe zugrunde, ein mit geringer 
Rechenzeit durchfiihrbares Verfahren zur automatischen Erzeu- 
gung einer Einhullenden, welche die Oberflache eines Kon- 
struktionsmodells mit vorgegebener Maximal abweichung approxi- 
miert, bereitzustellen. 

Die Aufgabe wird durch das Verfahren nach Anspruch 1 gelost . - 
Vorteilhafte Ausgestaltungen sind in den Unteranspruchen an- 
gegeben . 

Vorgegeben ist eine Vernetzung des Konstruktionsmodells . Die 
Vernetzung besteht aus Finiten Elementen mit Knotenpunkten . 
Das Verfahren umfafit folgende Schritte : 

- Ermittelt wird ein Quader, in dem die Vernetzung des Kon- 
struktionsmodells vollstandig enthalten ist. 

- Dieser Quader wird in Volumenelemente zerlegt. Die Zerle- 
gung ist eine vollstandige , so daS jeder Punkt des Quaders 
und damit jeder Punkt des Konstruktionsmodells in ein Vo- 
lumenelement f&llt. Der Quader wird so zerlegt, dafi jede 
Kante jedes dieser Volumenelemente kiirzer als oder genauso 
lang wie eine vorgegebene Schranke ist. Diese Schranke ist 
kleiner oder gleich der einzuhaltenden Maximal abweichung 
und kleiner oder gleich der fiir die nachf olgende Simulati- 
on geforderte Feinheit der Einhullenden. 

- Fiir jedes Volumenelement dieses Quaders wird gepruft, ob 
das Volumenelement mit mindestens einem ermittelten Fini- 
ten Element uberlappt . Eine Uberlappung liegt bereits vor, 
wenn das Finite Element und das Volumenelement einen ein- 
zigen Punkt gemeinsam haben. 
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. Die Menge der uberlappenden Volumenelemente bildet einen 
geometrischen Korper. Dieser Korper wird ermittelt. 

. Der Korper wird nach auSen durch die aufceren Begrenzungs- 
flachen der aufcenliegenden Volumenelemente begrenzt . Diese 
Begrenzungsflachen werden ermittelt. 

Die gesuchte Einhiillende wird aus der Menge der dergestalt 
ermittelten Begrenzungsflachen zusammengeset zt . 
Die erfindungsgemaS erzeugte Einhiillende weicht von der zu 
approximierenden Oberflache hochstens urn eine vorgegebene Ma- 
ximal abweichung ab. Die Zerlegung des Quaders, der das Kon- 
struktionsmodell vollstandig enthalt, wird mit einer vorgege- 
benen Genauigkeit durchgef uhrt . Die Zerlegung wird so durch- 
gefiihrt, da£ die groSte Kantenlange der Volumenelemente klei- 
ner oder gleich der Maximalabweichung ist. Diese groSte Kan- 
tenlange ist eine obere Grenze fur die Abweichung zwischen 
der Einhiillenden, die aus den auSenliegenden Begrenzungsfla- 
chen der uberlappenden Volumenelemente zusammengeset zt wird, 
und der durch die Vernet zung approximierten tatsachlichen 0- 
berflache des Konstruktionsmodells . Weil jede Kante kleiner 
oder gleich der vorgegebenen Schranke ist, ist die Maximalab- 
weichung hochstens so grofc wie diese vorgegebene Schranke. 
Denn die Kanten der Begrenzungsflachen sind hochstens so lang 
wie die maximale Kantenlange der Volumenelemente. Weil uber- 
lappende Volumenelemente ausgewahlt werden und keine anderen, 
stellt das Verfahren sicher, daS der Abstand zwischen der ap- 
proximierenden Vernetzung und der Einhiillenden hochstens die 
maximale Kantenlange und damit kleiner oder gleich der vorge- 
gebenen Maximalabweichung ist. 

Das Verfahren macht es iiberflussig, die vorgegebene Vernet- 
zung des Konstruktionsmodells zu verandern. Dadurch werden 
Rechenschritte und somit Verarbeitungszeit eingespart . Diese 
Einsparung ist urn so groSer, je f einer die vorgegebene Ver- 
netzung im Vergleich zur Maximalabweichung ist. Das Verfahren 
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erfordert bei N Quadern und M Finiten Elementen maximal M * N 
Priif schritte. 

Das erfindungsgemaSe Verfahren hangt nicht von der Art ab, 
wie das Konstruktionsmodell vernetzt wurde. Das Verfahren 
lafit sich insbesondere sowohl dann anwenden, wenn die Vernet- 
Z ung des Konstruktionsmodells aus Flachenelementen besteht, 
als auch dann, wenn sie aus Volumenelementen besteht. Das 
Verfahren lftSt sich auch fur jede Form der Flachenelemente 
anwenden. Der Quader wird vollig unabhangig von der Vernet- 
zung des Konstruktionsmodells in Volumenelementen zerlegt . 
Dadurch wird die Einhullende so erzeugt, wie die nachfolgende 
Simulation dies erfordert. Diese nachfolgende Simulation kann 
eine andere sein als diejenige, fur welche die vorgegebene 
Vernetzung erzeugt wurde . Beispielsweise wurde die Vernetzung 
fur eine Simulation der im System auftretenden Spannungen er- 
zeugt, die Einhullende wird hingegen fur eine Akustik- 
Simulation verwendet . 

Nur diejenigen Begrenzungsf lachen der uberlappenden Volumen- 
elemente werden ausgewahlt, die den durch die uberlappenden 
Volumenelemente gebildeten Korper nach auiSen begrenzen. Da- 
durch wird erreicht, da£ die innere Struktur des Konstrukti- 
onsmodells bei der Erzeugung der Einhullenden nicht beruck- 
sichtigt wird. Weiterhin wird erreicht, daS die erf indungsge- 
maJS erzeugte Einhullende ist daher immer geschlossen und re- 
gelmaKig ist . 

Falls das vorgegebene Konstruktionsmodell Locher aufweist, z. 
B. Bohrungen oder Taschen, so hangt es von deren Abmessungen 
ab, ob Approximationen dieser Locher in der Einhullenden auf- 
treten oder ob die Locher bei der Approximation nicht beruck- 
sichtigt werden. Falls der Durchmesser des Lochs groiSer ist 
als die groSte Kantenlange der Volumenelemente des Quaders, 
so liegt mindestens ein Volumenelement des Quaders ganz oder 
teilweise in dem Loch, und Begrenzungsf lachen dieses Quaders 
sind auSenliegende und gehoren zur Einhullenden. Ansonsten 
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wird das Loch vernachlassigt . Durch die frei wahlbare GroSe 
der Volumenelemente des Quaders wird also gesteuert, welche 
Locher beriicksichtigt und welche vernachlassigt werden. 
Die erzeugte Einhullende besitzt wesentlich weniger Knoten- 
punkte als die vorgegebene Vernetzung des Konstruktionsmo- 
dells. Eine Finite-Elemente-Simulation mit den Knotenpunkten 
der Einhullenden verwendet ein Gleichungssystem mit weniger 
Unbekannten und erf ordert daher weniger Rechenzeit . Diese 
Einsparung fallt insbesondere dann ins Gewicht, wenn die Sx- 
raulation mehrmals durchgefuhrt wird, z . B . urn verschiedene 
Konstruktionsstande miteinander zu vergleichen. 
Die Ausgestaltung nach Anspruch 2 erfordert besonders wenig 
Rechenzeit. Der Quader wird so zerlegt, daE die kvirzeste Kan- 
te jedes Volumenelements groSer oder gleich der langsten Kan- 
te jedes Finiten Elements der vorgegebenen Vernetzung ist. 
Die Zerlegung des Quaders ist frei wahlbar, und die Genauig- 
keit, mit der die Einhullende die Oberflache des Konstrukti- 
onsmodells approximieren soil, wird durch die Anf orderungen 
der Simulation, fur welche die Einhullende verwendet wird, 
vorgegeben. Daher ist es meistens moglich, daS die Volumen- 
elemente des Quaders langere Kanten als die Flachenelemente 
der Vernetzung des Konstruktionsmodells haben. 
Die Prufung, ob ein Finites Element der Vernetzung und ein 
Volumenelement des Quaders uberlappen oder nicht, wird durch 
die Ausgestaltung nach Anspruch 2 auf besonders einfache Wei- 
se durchgefuhrt. Gepruft wird, ob mindestens ein Knotenpunkt 
des Finiten Elements im Volumenelement liegt. Wird ein sol- 
dier Knotenpunkt gefunden, so wird die Prufung abgebrochen, 
und entschieden wird, daS das Volumenelement mit dem Finiten 
Element uberlappt. Falls kein Knotenpunkt des Finiten Ele- 
ments im volumenelement liegt, so kann auch kein anderer 
Punkt im Volumenelement liegen. Dies folgt aus der Ausgestal 
tung, dafi die Kanten des Volumenelements langer als die des 
Finiten Elements sind. Weil kein anderer Punkt des Finiten E 
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lements im Volumenelement liegt, uberlappen Finites Element 
und volumenelement nicht . Weil nur die Knotenpunkte unter- 
sucht werden und keine anderen Punkte, wird besonders wenxg 
Rechenzeit benotigt. 

im Folgenden wird ein Ausfuhrungsbei spiel der Erf indung naher 
beschrieben, dabei zeigen: 

Fig. 1. Knotenpunkte einer vorgegebenen unregelmaSigen Fla- 
che ; 

Fig. 2. ein die Flache enthaltendes Rechteck; 

Fig. 3. eine Zerlegung des Rechtecks in kleine Rechtecke; 

Fig. 4. die Rechtecke der Zerlegung, die Knotenpunkte enthal - 
ten; 

Fig. 5. die entstehende massive zweidimensionale Flache; 
Fig. 6. die erzeugte Einhullende. 

Das erfindungsgem&Se Verfahren wird zunachst vereinfacht an 
einem zweidimensionalen Beispiel veranschaulicht . In diesem 
Beispiel ist eine Einhullende fur eine unregelmalSige Flache 
in der Ebene zu erzeugen. Diese unregelmaSige Flache tritt in 
dieser Veranschaulichung an die Stelle des vorgegebenen Kon- 
struktionsmodells. Die zu erzeugende Einhullende ist ein 
Streckenzug in dieser Ebene. Fig. 1 zeigt einige der Knoten- 
punkte derjenigen Flachenelemente, die zu einer Vernetzung 
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der Flache gehoren. Der ubersichtlichkeit halber sind viele 
Knotenpunkte nicht dargestellt. Die durch die Knotenpunkte 
gebildeten Flachenelemente der unregelmaSigen Flache sind e- 
benfalls nicht dargestellt. 

Bin Rechteck wird ermittelt, in dem die vorgegebene unregel- 
mafiige Flache und damit die Vemetzung enthalten ist. Die un- 
regelmaSige Flache und das ermittelte Rechteck werden in 
Fig. 2 gezeigt. Dieses Rechteck tritt an die Stella des Qua- 
ders . 

Das Rechteck wird in gleich groSe Rechtecke mit vorgegebenen 
Kantenlangen zerlegt. Fig. 3 zeigt diese Rechtecke der Zerle- 
gung, die an die Stelle der Volumenelemente treten, sowie die 
Knotenpunkte der Fig. 1. 

Gepriift wird, welche Rechtecke der Zerlegung jeweils mindes- 
tens einen Knotenpunkt der unregelmaSigen Flache umfassen. 
Diese iiberlappenden Rechtecke sind in Fig. 4 schraffiert dar- 
gestellt . 

Die in Fig. 4 schraffiert dargestellten Rechtecke bilden eine 
zweidimensionale Flache. Die Hohlraume in dieser zweidimensi- 
onalen Flache werden ausgef ullt . Dadurch entsteht die in 
Fig. 5 gezeigte massive zweidimensionale Flache. 
Die Begrenzungslinien dieser massiven Flache werden zur ge- 
suchten Einhullenden zusammengeset zt . Diese Begrenzungsli- 
nien, die einen Streckenzug bilden, werden in Fig. 6 durch 
fette Linien gezeigt. 

Das nun folgende Ausf iihrungsbei spiel bezieht sich auf einen 
Bestandteil eines neuen Kraf tf ahrzeuges als dem technischen 
System. Die erzeugte Einhullende wird beispielsweise dazu 
verwendet, die Schallemission von einem festen Korper, z. B. 
einem Motor Oder einer Karosserie eines Kraf tf ahrzeuges , 
durch eine Finite-Elemente-Simulation zu untersuchen. Die 
Schallemission wird durch Korperschall hervorgeruf en. Dieser 
entsteht durch Zug- und/oder Schubspannungen in Festkorpern. 
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Die Schall emission wird durch Vibrationen an der Oberflache 
des Motors oder der Karosserie hervorgeruf en . Diese Vibratio- 
nen wirken auf die umgebende Luft ein. Die Vibrationen werden 
in einen Anteil senkrecht zur Oberflache und einen parallel 
zur Oberflache zerlegt . Nur die Anteile der Vibrationen senk- 
recht zur Oberflache werden in der Akustik-Simulation beruck- 
sichtigt, nicht aber die Anteile parallel zur Oberflache. 

Aus der Strukturdynamik des Motors und seiner Oberflache re- 
sult iert eine Verteilung der Schnelle der Vibrationen. Diese 
Schnelle wird durch eine Akustik-Simulation ermittelt. Mit 
der Schnelle wird die Auslenkungsgeschwindigkeit ds/dt be- 
zeichnet, wobei s(t) die Auslenkung der Schallwelle senkrecht 
zur Oberflache in Abhangigkeit von der Zeit t ist und wobei 
gilt: 

ds/dt * -co * A * sin (ot + kx + <j>) . 

Die Strukturdynamik des Motors wird mit einer ersten Finite- 
Elemente-Simulation berechnet . Die berechnete Strukturdynamik 
umfaSt die Verschiebung, die Geschwindigkeit oder die Be- 
schleunigung z. B. der Verschiebung der Oberflache. Aus die- 
ser wird die Schallschnelle der Oberflache berechnet, 

Vorgegeben ist ein rechnerverf ugbares Konstruktionsmodell ei- 
nes zu konstruierenden oder zu untersuchenden Systems, z. B. 
des Motors. Dieses Konstruktionsmodell ist auf einem Daten- 
speicher einer Datenverarbeitungsanlage abgespeichert und hat 
vorzugsweise die Form eines dreidimensionalen CAD-Modells. 
Das Konstruktionsmodell wird beispielsweise mit Hilfe von 
Punkten, Kurven, Vektoren, Linien, Konturen, polygoniale Net- 
ze, gekrummten Oberflachen und/oder Volumenelementen in dem 
Datenspeicher reprasentiert . Das erf indungsgemafie Verfahren 
wird mit dieser Datenverarbeitungsanlage oder einer anderen, 
die Lesezugriff auf den Datenspeicher hat, durchgef uhrt . 

Die Methode der Finiten Elemente ist aus „Dubbel - Taschen- 
buch fur den Maschinenbau" , 20. Auflage, Springer- Verlag, 
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2001, C 48 bis C 50, aus B. Klein: „FEM - Grundlagen und An- 
wendungen der Finite-Elemente-Methode" , Vieweg-Verlag , 
3. Auflage, 1999, aus T. R. Chandrupalta & A. D. Belegundu: 
„ Introduction to Finite Element in Engineering", Prentice- 
Hall, 1991, sowie aus DE 19927941 CI bekannt . Durch eine Fi~ 
nite-Elemente-Simulation werden physikalische Vorgange unter- 
schiedlicher Arten untersucht, beispielsweise Schwingungs vor- 
gange, z. B . die Gerauschemission von der Oberflache eines 
Motors oder einer Karosserie eines Kraf tf ahrzeuges . 

In dem vorgegebenen Konstruktionsmodell wird eine bestimmte 
Menge von Punkten festgelegt, die Knotenpunkte heiSen* Als 
Finite Elemente werden diejenigen Flachen- oder Volumenele- 
mente bezeichnet, deren Ecken durch Knotenpunkte definiert 
werden. Die Knotenpunkte bilden ein Netz in dem Konstrukti- 
onsmodell, weswegen der Vorgang, Knotenpunkte festzulegen und 
Finite Elemente zu erzeugen, Vernetzen und die Finiten Ele- 
mente Vernetzung des Modells genannt werden. In DE 
10010408 Al und DE 19 933314 Al werden Verfahren zum Vernetzen 
eines gegebenen Konstruktionsmodells beschrieben . 

Moglich ist, das Konstruktionsmodell beim Vernetzen vollstan- 
dig in Finite Elemente zu zerlegen. In diesem Fall gehort je- 
der Punkt des Konstruktionsmodells zu mindestens einem Fini- 
ten Element. Oft reicht aber eine Naherungslosung aus, die 
Knotenpunkte einspart und damit zu einem kleineren Glei- 
chungssystem f uhrt . Falls beispielsweise die Karosserie eines 
Kraf tf ahrzeugs untersucht wird, so reicht es meist aus, die 
Blechteile der Karosserie durch Flachenelemente naherungswei- 
se zu vernetzen. Fugeverbindungen zwischen Blechteilen, z. B. 
SchweiE- oder Klebenahte, werden mittels Volumenelementen 
vernetzt. Dadurch lassen sich Vorgange in den Fugeverbindun- 
gen, z. B. Verschiebungen und Querspannungen, genauer unter- 
suchen. In manchen Anwendungen reicht es aus, die Oberflache 
des Konstruktionsmodells zu vernetzen, beispielsweise durch 
eine „Triangulierung w mittels dreieckiger Flachenelemente. 
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Das erf indungsgemafie Verfahren laSt sich ohne Modifikation 
fur jede dieser Vernetzungen anwenden. Eine beliebige Vernet- 
zung des Konstruktionsmodells wird dem Verfahren vorgegeben. 

Je nach Pr obi ems t el lung werden die Verschiebungen dieser Kno- 
tenpunkte und/oder Rotationen der Finiten Elemente in diesen 
Knotenpunkten und/oder die Spannungen in bestimmten Punkten 
dieser Finiten Elemente, namlich in den Integrationspunkten, 
als Unbekannte eingefiihrt. Gleichungen werden aufgestellt, 
welche z. B. die Verschiebungen, Rotationen und/oder Spannun- 
gen innerhalb eines Finiten Elements naherungsweise beschrei- 
ben. Weitere Gleichungen result ieren aus Abhangigkeiten zwi- 
schen verschiedenen Finiten Elementen. Insgesamt wird ein oft 
sehr umf angreiches Gleichungssystem mit den Knotenpunkt- 
Verschiebungen, Knotenpunkt -Rotationen, Spannungen in Integ- 
rationspunkten oder weitere Grofeen als Unbekannte aufgestellt 
und numerisch gelost. Aus der Losung des Gleichungs systems 
lassen sich Aussagen iiber das physikalische Verhalten des 
Systems gewinnen. 

Die Einhullende laSt sich auch fur eine Simulation gemaS der 
Boundary-Elemente-Methode verwenden. Diese Methode wird in 
Dubbel, a.a.O., C50 - C52, beschrieben. Wie eine Akustik- 
Simulation mit der Finite-Elemente-Methode oder der Boundary- 
El emente- Methode durchgefuhrt wird, ist z. B. in Dubbel, 
a.a.O., 033 - 035, beschrieben. 

Zur Durchfuhrung des erf indungsgemaSen Verfahrens wird ledig- 
lich die Vernetzung des Konstruktionsmodells benotigt . Weite- 
re Informationen iiber das Konstruktionsmodell , z. B. Kriim- 
mungsradien oder Fasen oder die Oberf lachenbeschaf f enheit , 
erfordert das Verfahren nicht. Die Datenverarbeitung, die das 
Verfahren durchfuhrt, liest die Vernetzung des Konstruktions- 
modells ein. Hierbei werden die im Folgenden beschriebenen 
Informationen iiber Knotenpunkte und iiber die Finiten Elemente 
eingelesen. 
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Ein Knotenpunkt ist durch eine eindeutige Kennung, z. B. eine 
fortlaufende Nummer, und seiner Position gekennzeich.net . Die 
Kennung ist in der Vernetzung eindeutig, weswegen zwei ver- 
schiedene Knotenpunkte stets zwei verschiedene Kennungen ha- 
ben. Die Position des Knotenpunkts wird durch je eine x- , y- 
und z-Koordinate in einem dreidimensionalen Koordinatensystem 
gekennzeichnet . Beispielsweise um Verschiebungen und Verfor- 
mungen unter dem Einflufi eines Lastfalls zu untersuchen, kann 
vorgesehen sein, dafi zwei Knotenpunkte zweier verschiedener 
Finiten Elemente vor der Verformung dieselbe Position und - 
daher dieselben x~ , y- und z-Koordinaten haben. Diese beiden 
Knotenpunkte werden durch verschiedene Kennungen unterschie- 
den. 

Man unterscheidet zweidimensionale Finite Elemente, die im 
folgenden Flachenelemente genannt werden, und dreidimensiona- 
le Finite Elemente, die im Folgenden Volumenelemente heifeen. 
Gekrummte Flachen und Korper, die durch gekrummte Flachen im 
Raum angenahert werden, werden oft mittels Flachenelementen 
vernetzt. Die gebrauchlichsten Flachenelemente sind Dreiecke 
und Vierecke. Alle Ecken der Flachenelemente sind Knotenpunk- 
te. Manchmal werden Vierecke mit sechs Knotenpunkten verwen- 
det . Neben rechteckigen Flachenelementen konnen auch andere 
viereckige Flachenelemente vorgesehen sein. Flachenelemente 
werden vorzugsweise durch eine isoparametrische Darstellung 
beschrieben, die z. B. aus T. R » Chandrupalta & A. D. Bele- 
gundu, a.a.O., bekannt ist. 

Als Volumenelemente werden in der Regel Wurfel und andere 
Quader verwendet, also Volumenelemente mit sechs ausschlieS- 
lich rechteckigen Begrenzungsf lachen. Derartige Volumenele- 
mente werden oft als Hexaeder bezeichnet . Typischerweise sind 
die acht Ecken des Quaders Knotenpunkte des Volumenelements . 
Verwendet werden Volumenelemente mit acht bis zwanzig Knoten- 
punkten. 
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Jedes Finite Element ist durch eine in der Vernetzung eindeu- 
tige Kennung, z. B. eine fortlaufende Nummer, und durch seme 
Knotenpunkte gekennzeichnet . Vorzugsweise umfassen die infor- 
mations uber das Finite Element eine Auflistung der Kennun- 
gen seiner Knotenpunkte sowie eine Kennzeichnung, von welcher 
Art das Finite Element ist. 

Als ersten Schritt bei der Durchfiihrung des Ausfuhrungsbei- 
spiels liest die Datenverarbeitung, die das Verfahren aus- 
fuhrt, die Vernetzung ein. Ein Quader wird ermittelt, in dem 
die Vernetzung des Konstruktionsmodells vollstandig enthalten 
ist. Vorzugsweise wird diese Quader so gelegt, daS jede sei- 
ner Begrenzungsflachen senkrecht auf einer Koordinatenachse 
liegt und somit seine Kanten achsenparallel ausgerichtet 
sind. Jede der sechs Begrenzungsflachen liegt dann in einer 
der folgenden Ebenen: 

{ (x,y,z) | x = x_min }, { (x,y,z) | x = x_max }, 
{ (x,y,z) | y = y_min }, { (x,y,z) | y = y_max }, 
{ (x,y,z) | z = z_min }, { (x,y,z) | z = z_max }. 
Der Wert x_min wird mit Hilfe der Positionen der Knotenpunkte 
der eingelesenen Vernetzung berechnet, namlich als der 
kleinste Wert aller x-Koordinaten der Knotenpunkte. Entspre- 
chend wird x_max als der groSte Wert aller x-Koordinaten der 
Knotenpunkte berechnet . Auf analoge Weise werden y_min, 
y_max, z_min und z_max berechnet. Alle Finiten Elemente der 
vorgegebenen Vernetzung liegen innerhalb dieses Quaders, well 
jeder Eckpunkt eines Finiten Elements ein Knotenpunkt der 
Vernetzung ist und weil die Berechnung von x_min, . . . , z_max 
sicherstellt, da£ alle Knotenpunkte innerhalb des Quaders 
liegen. 

Der Quader, in dem die Vernetzung enthalten ist, wird in Vo- 
lumenelemente zerlegt. Vorzugsweise haben alle Volumenelemen- 
te die gleichen Abmessungen und die gleiche geometrische Form 
und grenzen nahtlos aneinander. Beispielsweise wird der Qua- 
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der in gleich grofee Wiirfel zerlegt. In dem Falle, da£ ein 
harmonischer Schwingungsvorgang durch eine Finite-Elemente- 
Simulation auf Basis der Einhullenden untersucht werden soli, 
wird wie folgt festgelegt, wie der Quader in Wiirfel zerlegt 
wird: Eine halbe Welle des Schwingungsvorgangs wird in der 
Finite-Elemente-Simulation durch einen Streckenzug mit drei 
gleich langen Strecken approximiert , eine voile Welle demnach 
durch sechs Strecken- Daher wird der Quader in Wiirfel zer- 
legt, deren Kantenlange ein Sechstel der Wellenlange, also 
der Lange der vollen Welle, ist, Diese Wellenlange hangt von 
der Frequenz des Schwingungsvorgangs und von dem Medium, in 
dem die Schwingung sich ausbreitet, ab, und ist daher be- 
kannt . 

Gepruft wird, ob die kiirzeste Kantenlange der Volumenelemente 
des Quaders grower ist als die langste Kante der Vernetzung 
des Konstruktionsmodells . Ist dies nicht der Fall, so wird 
die Vernetzung des Konstruktionsmodells verf einert . 

Vorzugsweise werden zunachst die Kantenlangen der Volumenele- 
mente in x-, y- und z-Richtung ermittelt oder vorgegeben. An- 
schlieSend werden die Langen der drei Kanten des Quaders in 
x- , y- und z-Richtung durch die vorgegebenen oder ermittelten 
Kantenlangen der Volumenelemente dividiert . Die drei Quotien- 
ten werden auf die nachsten natiirlichen Zahlen N_x, N__y und 
N_z auf gerundet . Dadurch wird die Anzahl der Volumenelemente 
in x-, y- und z-Richtung berechnet . Eine tatsachliche Kanten- 
lange eines Volumenelements kann daher von der vorgegebenen 
leicht abweichen. 

Die Volumenelemente, in die der Quader zerlegt wird, werden 
fortlaufend durchnumeriert . Dadurch erhalt jedes dieser Volu- 
menelemente eine eindeutige Kennung. Vorzugsweise ist die 
Kennung ein Tripel (n_x, n_jy, n_z) , das die relative Position 
des Volumenelements im Quader angibt . Das Volumenelement , das 
den Eckpunkt (xjnin, y_min, z_min) des Quaders enthalt, er- 
halt die Kennung (1, 1, 1) . Das Volumenelement, das den Eck- 
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punkt (x_max, y_max , z_max) des Quaders enthalt, erhalt die 
Kennung (N_x, N_y, N_z) , wobei N_x, N__y und N_z die Anzahl 
der Volumenelemente in x- , y~ bzw. z-Richtung angibt . 

Als nachstes wird ermittelt, welche Volumenelemente des Qua- 
ders mit Finiten Elementen der Vernetzung uberlappen. Hierzu 
wird eine Datenstruktur in Form einer dreispaltigen Tabelle 
erzeugt . Anstelle einer Tabelle kann auch eine verkettete 
Liste als Datenstruktur verwendet werden, wobei jedes Listen- 
element drei Eintrage ermoglicht. 

In der ersten Spalte der Tabelle werden die Kennungen der Vo- 
lumenelemente des Quaders eingetragen, in die zweite Spalte 
in jeder Zeile eine Null. AnschlieSend wird eine Datenstruk- 
tur, in der die Knotenpunkte der Vernetzung abgespeichert 
sind, durchlaufen. Fur jeden Knotenpunkt der Vernetzung wird 
ermittelt, in welches Volumenelement des Quaders dieser Kno- 
tenpunkt fullt. Wie oben dargelegt, liegt jeder Knotenpunkt 
in dem Quader, und der Quader wird vollstandig und ohne Rest 
in Volumenelemente zerlegt . Moglich ist, da£ ein Knotenpunkt 
auf der Begrenzungsf lache zweier aneinander angrenzenden Vo- 
lumenelement en liegt. In der zweispaltigen Tabelle wird fur 
die Volumenelemente, in die der Knotenpunkt fallt, der alte 
Eintrag in der zweiten Spalte um Eins vergroiSert . Nach Ende ~ 
des Schrittes gibt die Tabelle fur jedes Volumenelement an, 
wie viele Knotenpunkte in ihm liegen. 

Die Menge der uberlappenden Volumenelemente ist genau die 
Menge derjenigen Volumenelemente, die mindestens einen Kno- 
tenpunkt enthalten und in der zweiten Spalte der Tabelle ei- 
nen Wert groSer Null besitzen. Diese Menge bildet einen drei- 
dimensionalen Korper, der vollstandig im Quader liegt. Dieser 
Korper kann Hohlraume enthalten, die im nachsten Schritt ge- 
fullt werden. Nach Auf f til lung der Hohlraume entsteht ein mas - 
siver Korper. Die dritte Spalte der Tabelle legt fest, welche 
Volumenelemente des Quaders zu diesem massiven Korper gehoren 
und welche nicht- Hierzu wird in die dritte Spalte der Tabel- 
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le eine Eins eingetragen, wenn in der zweiten Spalte ein Wert 
groSer Null steht . 

Die Volumenelemente, die keinen Knotenpunkt enthalten und da- 
her in der zweiten Spalte eine Null haben, gehoren entweder 
zu Hohlraumen im Korper gehoren oder liegen aufeerhalb des 
Korpers . Um diese beiden Falle zu unterscheiden, wird wie 
folgt vorgegangen: Fur jedes Volumenelement VE, das keinen 
Knotenpunkt enthalt und fiir das in der dritten Spalte noch 
kein Eintrag vorgenommen wurde, werden die sechs angrenzenden 
Volumenelemente ermittelt, also diejenigen, die mit einer Be- 
grenzungsf lache an VE angrenzen. Ermittelt wird, welche die- 
ser sechs angrenzenden Volumenelemente mindestens einen Kno- 
tenpunkt enthalten und welche keinen. Fur die an VE angren- 
zenden Volumenelemente, die keinen Knotenpunkt enthalten, 
werden wiederum die angrenzenden Volumenelemente ohne Knoten- 
punkt e ermittelt, die nicht an VE angrenzen. Diese Ermittlun- 
gen werden fortgesetzt, bis entweder ein zusammenhangender 
Hohlraum, der VE enthalt, ermittelt wurde oder bis der Rand 
des Quaders erreicht ist . Alle Volumenelemente eines zusam- 
menhangenden Hohlraums werden dadurch „gefullt", daiS in der 
dritten Spalte jeweils eine Eins eingetragen wird. Wird hin- 
gegen ausgehend von VE der Rand des Quaders erreicht, so 
liegt VE nicht im Inneren des Korpers, und in der dritten 
Spalte werden fur VE und fur die angrenzenden Volumenelemente 
ohne Knotenpunkte je eine Null eingetragen. 

Das Auffiillen von Hohlraumen laiSt sich auch auf anderem Wege 
durchfuhren. Zunachst wird in der dritten Spalte fur alle Vo- 
lumenelemente mit je mindestens einem Knotenpunkt je eine 
Eins eingetragen. Anschliefeend werden nacheinander alle Volu- 
menelemente ermittelt, die am Rande des Quaders liegen, d. h. 
die mindestens eine Begrenzungsf lache besitzen, die in einer, 
Begrenzungsf lache des Quaders liegt. Diejenigen dieser Volu- 
menelemente, die keinen Knotenpunkt der Vernetzung enthalten, 
sind au£enliegende Volumenelemente und erhalten eine Null in 
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der dritten Spalte. Diejenigen Volumenelemente, die an diese 
Volumenelemente am Rande angrenzen und ebenfalls keine Kno- 
tenpunkte enthalten, werden ermittelt und erhalten ebenfalls 
eine Null in der dritten Spalte. Wiederholt werden angrenzen- 
de Volumenelemente ohne Knotenpunkte ermittelt, bis alle Vo- 
lumenelemente untersucht sind und die dritte Spalte vollstan- 
dig ausgefullt ist. 

Nach dem Auffullen liegt eine aus Volumenelementen des Qua- 
ders bestehende Vernetzung eines massiven dreidimensionalen 
Korpers vor. Diejenigen Begrenzungsf lachen dieser Vernetzung 
werden ermittelt, die diesen Korper nach auSen begrenzen. 
Vorzugsweise wird hierzu aus dem Volumenmodell , das ist die 
Vernetzung des massiven Korpers mit Volumenelementen, ein O- 
berf lachenmodell , das ist eine Vernetzung der Oberf lache mit 
Flachenelementen, erzeugt . Dieser Schritt wird z. B. mit ei- 
nem Praprozessor fur Finite-Elemente-Werkzeuge durchgef uhrt . 
Als erstes wird eine Beschreibung der aus Volumenelementen 
bestehenden Vernetzung des Korpers in einer Datei abgespei- 
chert, z, B. im Datenformat „NASTRAN Bulk Data". Der Prapro- 
zessor „MEDINA U vermag eine Vernetzung aus einer Datei einzu- 
lesen, wenn diese in einem standardisierten Datenformat , z. 
B. in NASTRAN Bulk Data, vorliegt. Eine Beschreibung von 
MEDINA ist unter 

http : //www. c3pdm. com/des/products/medina/documentation/medina 
-DINA4_e.pdf, abgefragt am 27. 11. 2003, verfugbar, MEDINA 
besitzt eine Funktionalitat , urn aus einem Volumenmodell ein 
Oberf lachenmodell zu erzeugen. Dieses Oberf lachenmodell ist 
die gesuchte Einhullende. 

Vorzugsweise wird das Oberf lachenmodell geglattet („rela- 
xiert") , wobei Treppenstuf en in der Vernetzung entfernt wer- 
den. Auch die Glattung wird z. B. mit einer entsprechenden 
Funktionalitat von MEDINA durchgef uhrt . 

Die erf indungsgemafi erzeugte Einhullende laSt sich auch dazu 
verwenden, Ergebnisse der Finite-Elemente-Simulation, z. B. 
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das Schwingungsverhalten oder die Schallausbreitung, zu visu- 
alisieren. Fur diese Visual! si erung werden nur die berechne- 
ten physikalischen GroiSen an der Oberflache des untersuchten 
Systems benotigt, nicht aber die innere Struktur. Mit der 
Einhullenden wird eine Datenreduktion durchgef uhrt . Die Be- 
rechnungsergebnisse werden auf die Flachenelemente der Ein- 
hiillenden umgerechnet, z. B. indem per Interpolation die Wer- 
te der Grofeen in den Eckpunkten dieser Flachenelemente be- 
rechnet werden. Die berechneten Grofien werden vorzugsweise 
auf einem Bildschirm einer Datenverarbeitungsanlage oder auf 
einem Papierausdruck visualisiert . Infolge der Datenreduktion 
lafit sich die Visual isierung schneller erzeugen und veran- 
dern, z. B. wenn ein Benutzer die Visual isierung auf dem 
Bildschirm drehen oder vergroSern will. 
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19. 12 .2003 



Patent anspruche 



1. Verfahren zur automatischen Erzeugung einer Einhullenden, 
welche die Oberflache eines reclmerverf ugbaren Konstruk- 
tionsmodells eines technischen Systems approximiert , 

wobei eine aus Finiten Elementen mit Knotenpunkten beste- 
hende Vernetzung des Konstruktionsmodells vorgegeben ist 

und das Verfahren die Schritte umf a£t : 

- Ermitteln eines Quaders, in dem die Vernetzung voll- 
standig enthalten ist, 

Zerlegen des Quaders in Volumenelemente mit Kanten, 
die hochstens genauso lang wie eine vorgegebene 
Schranke sind, 

Fur jedes Volumenelement Priifen, ob das Volumenelement 
mit mindestens einem Finiten Element der Vernetzung u- 
berlappt , 

- Ermitteln eines Korpers, der durch die liberlappenden 
- Volumenelemente gebildet wird, 

- Ermitteln derjenigen Begrenzungsf lachen, die den Kor- 
per nach auEen begrenzen, und 
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Zusammensetzen der Einhullenden aus den ermittelten 
Begrenzungsf lachen . 

2 . Verf ahren nach Anspruch 1 , 
dadurch gekennzeichnet , 
da£ 

der Quader so zerlegt wird, da£ die kurzeste Kante je~ 
des Volumenelements grower oder gleich der langsten 
Kante jedes Finiten Elements ist 

dann, wenn es einen Knotenpunkt eines Finiten Elements 
gibt, der in einem Volumenelement liegt, Entscheiden, 
daS das Volumenelement mit dem Finiten Element uber- 
lappt, und 

dann, wenn kein Knotenpunkt eines Finiten Elements in 
einem Volumenelement liegt, Entscheiden, da£ dieses 
Volumenelement mit dem Finiten Element nicht uber- 
lappt . 

3 . Verf ahren nach Anspruch 1 oder Anspruch 2 , 
dadurch gekennzeichnet , daS 

die Finiten Element e der Vernetzung Rechtecke 
und die Volumenelemente Quader sind. 

4. Verf ahren nach einem der Anspriiche 1 bis 3, 
dadurch gekennzeichnet, 
da£ 

eine Wellenlange eines harmonischen Schwingungsvor- 
gangs vorgegeben ist, 
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als Volumenelemente Wurfel verwendet werden, deren 
Kantenlange ein vorgegebener Bruchteil der Wellenlange 
ist, und 

unter Verwendung der Einhullenden eine Finite- 
Element e- Simulation des Schwingungsvorgangs durchge- 
fiihrt wird. 

5. Verfahren nach einem der Anspruche 1 bis 4, 
dadurch gekennzeichnet , 

da£ beim Ermitteln der Begrenzungsf lachen 

diejenigen Volumenelemente des Quaders ermittelt wer- 
den, die Hohlraume im Inneren des ermittelten Korpers 
bilden, 

diese Hohlraume gefiillt werden und 

die Begrenzungsf lachen des dadurch erzeugten massiven 
dreidimensionalen Korpers ermittelt werden. 

6. Computerprogramm-Produkt , das direkt in den internen 
Speicher eines Computers geladen werden kann und Soft- 
wareabschnitte umfaSt, mit denen ein Verfahren nach einem 
der Anspruche 1 bis 5 ausgefiihrt werden kann, wenn das 
Produkt auf einem Computer lauft. 

7. Computerprogramm-Produkt, das auf einem von einem Compu- 
ter lesbaren Medium gespeichert ist und das von einem 
Computer lesbare Programm-Mittel aufweist, die den Compu- 
ter veranlassen, ein Verfahren nach einem der Anspruche 1 
bis 5 auszufuhren. 
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Meyer -Gramann 

DaimlerChrysler AG 

19 . 12 .2003 



Zusammenf as sung 



Die Erfindung betrifft ein Verfahren zur automatischen Erzeu- 
gung einer Einhullenden fur ein rechnerverfugbares Konstruk- 
tionsmodell, das durch eine vorgegebene Vemetzung in Finite 
Elemente mit Knotenpunkten zerlegt ist. Ermittelt wird ein 
Quader, in dem die Vernetzung vollstandig enthalten ist. Die- 
ser Quader wird in Volumenelemente zerlegt . Die Zerlegung ist 
eine vollstandige , so daE jeder Punkt des Quaders und damit 
jeder Punkt des Konstruktionsmodells in ein Volumenelement 
fallt. Fur jedes Volumenelement dieses Quaders wird gepruft, 
ob das Volumenelement mit mindestens einem ermittelten Fini- 
ten Element iiberlappt . Eine Uberlappung liegt bereits vor, 
wenn das Finite Element und das Volumenelement einen einzigen 
Punkt gemeinsam haben. Die Menge der uberlappenden Volumen- 
elemente bildet einen geometrischen Korper. Dieser wird nach 
auiSen durch die auSeren Begrenzungsflachen der auSenliegenden 
Volumenelemente begrenzt . Diese Begrenzungsflachen werden er- 
mittelt. Die gesuchte Einhullende wird aus der Menge der der- 
gestalt ermittelten Begrenzungsflachen zusammengeset zt . Das 
Verfahren arbeitet schnell und mit vorgegebener Genauigkeit. 

(Fig. 4) 
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